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Abstract We review modeling attempts for the paradigmatic contact process 
(or SIS model) on adaptive networks. Elaborating on one particular proposed 
mechanism of topology change (rewiring) and its mean field analysis, we obtain 
a coarse-grained view of coevolving network topology in the stationary active 
phase of the system. Introducing an alternative framework applicable to a wide 
class of adaptive networks, active stationary states are detected, and an extended 
description of the resulting steady-state statistics is given for three different 
rewiring schemes. We find that slight modifications of the standard rewiring 
rule can result in either minuscule or drastic change of steady-state network 
topologies. 



1 Introduction 

In the past two decades, complex networks - multiple agents (nodes) engaging in structured 
interactions (links connecting nodes) - have become a prominent paradigm in tackling com- 
plex systems HI2I3I41 . More recently, allowing for the coevolution of network structure and 
node-state dynamics has given rise to abundant literature on adaptive networks 15I6L captur- 
ing such diverse phenomena as the emergence of cooperation [7 8 9 10 1 1I12L opinion for- 
mation 113I14I15I16I1 . disease spreading II17I18I19I20I21I . speciation 1221231 and traffic flows 
[24 25 26]. While some contributions explore the respective phenomenology with individual- 
based simulations [8 10 27 28), others also focus on providing explanatory frameworks for 
observed dynamics 071911 II 151 191201291 . 

Cyclic processes on adaptive networks, i.e. nodes going through a cyclic sequence of 
states while changing their connections, may yield adaptive networks displaying perpetual 
node-state and link dynamics Q. Those processes correspond to the active phase of the sys- 
tem, as opposed to the frozen phase where the system reaches a static equilibrium. As system 
dynamics often rely on the existence of active links connecting nodes of different states, the 
frozen phase either manifests itself in a globally polarized network (with all nodes carrying 
the same state) or a fragmented network (featuring two or more connected components that 
are polarized, whereas globally the network is not). Transitions to and between those ab- 
sorbing states in adaptive networks have been dealt with using mean-field approximations or 
small perturbations of the fragmented state 1151181301311321331 . 

Node-state and link dynamics in the active phase of adaptive networks can be highly com- 
plex 134 351, and stationary (as well as oscillatory) steady states have been shown to occur 
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for a variety of frameworks both in Monte-Carlo (MC) simulations and in their approximate 
mean field descriptions |6). For such networks in dynamic equilibrium (DE), nodes undergo 
permanent state and degree evolution, while global statistics characterizing ensemble dynam- 
ics and network topology settle down to a steady state. Describing network statistics in DE is 
essential to understanding the relation between dynamics and network structure and to apply 
these ideas to real world examples. 

Adaptive networks in the active phase have been treated analytically when timescales of 
node dynamics and topology evolution are separable 1 36 37 38]. For fast network dynamics 
in the context of evolutionary games, it was shown that the node states evolve according to an 
effective pay-off matrix that takes into account the equilibrium network properties 061371 . 

For the more general case that permits similar timescales, two equation-based frameworks 
taking the contact process as the underlying dynamics have so far been put forward. One is 
based on the pairwise formalism of |fl9l and related frameworks I20I21I39L where the time 
evolution of network motif densities is modeled up to the level of pairs, with description of 
highest-order motifs relying on the standard pair approximation moment closure assumption. 
This yields a low-dimensional set of ordinary differential equations (ODEs) that allows for 
analytic treatment, predicting system phases and pair densities in transients and DE with an 
accuracy limited by moment closure validity. Approximate phase diagrams were also derived 
in the scope of even simpler effective mean-field descriptions 1391 . While this type of models 
gives a global account of the network in DE through averages that characterize steady-state 
dynamics, its low number of degrees of freedom precludes any detailed description of the 
network's topology. The approach used in Ifl9ll for the contact process was extended to a 
three-state model (the Susceptible-Infected-Recovered-Susceptible cycle) describing more 
realistic infection dynamics l40l . Both models have stationary active phases, and simulations 
in DE converge to well-defined overall and state-specific degree distributions. In l40l . we 
also find an attempt to describe the topology of the network by translating the pairwise dy- 
namics of (T9| to a degree-class formulation similarly to the ansatz in l4lTl . This approach 
yields a self-contained method to determine the degree distributions of all three node types, 
but fails to reproduce the observed output of MC simulations. Given however the infected de- 
gree distribution extracted from simulations, the two remaining distributions are accurately 
described. 

Elaborating on a state- and degree-based compartmental formulation developed earlier 
l42l . Marceau et al. in 1431 on the other hand modeled the time evolution of fractions of 
nodes with the same status and joint degree (in a network with n possible node states, the 
joint degree (k\, fe, k n ) of a node is the set of numbers kj of its connections to nodes 
of type j G {1,2, n}). In the spirit of the moment closure in the aforementioned ansatz, 
infection dynamics beyond a node's immediate neighborhood are approximated by a mean 
field that is computed en route and assumes no correlations beyond the level of next neigh- 
bors. A large set of coupled nonlinear ODEs ensues that defies analytic treatment. Instead, 
numerical integration yields the time evolution of the joint degree distribution, as well as 
of low-order network motif densities derived from it. This approach provides an alternative 
to stochastic simulations on networks for the description of the DE, with the same limita- 
tions due to transient's length and the additional loss of accuracy involved in the mean field 
approximation. 

A third approach goes back to the original formulation of [19l and avoids any moment 
closure approximation in the numerical integration of the ODEs by computing the triplet 
densities at each integration step through short bursts of MC simulations on networks l44l . 
Much like 1401 . this hybrid approach gives a more precise ODE-based description of the 
system than its purely deterministic counterpart. It accurately reproduces the time evolution 
of the network's node and pair densities, and hence also the global phase diagram, but it 
brings no improvement regarding the analysis of the steady-state degree distributions. 

Individual-based simulations of various adaptive networks in DE show that apart from 
global averages and degree-related probability distributions, a wide range of other topologi- 
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cal measures settle down to an equilibrium. The most comprehensive description of DEs in 
frameworks like II43I45I is at the level of configuration model networks, implying that only 
the topology measures determined by the joint degree distributions can be derived. To this 
point a more extensive account of steady-state topology has to resort to modeling the time 
evolution of the network's adjacency matrix or a related construct 11461 . In 11471 however, the 
steady-state community structure of a general class of adaptive networks in DE was charac- 
terized by a simple ODE for the network's modularity. 

In this paper, we will focus on a popular adaptive system: the Susceptible-Infected- 
Susceptible (SIS) model with rewiring by Gross et al. (T9 |. For that particular rewiring sce- 
nario, more predictions from their pairwise model are extracted. A novel analytic framework 
that stochastically models adaptive networks in DE is then introduced and applied to the con- 
tact process with rewiring, elaborating on results presented in l45l . Finally, different rewiring 
rules are considered, which in SIS terminology represent disease awareness under different 
conditions, and the proposed framework is used to show how the degree distributions in DE 
depend on them. 



2 SIS with Selective Rewiring in the Pairwise Model 

As a simple cyclic process, the contact process on an adaptive network was proposed in |fl9l 
to model the spreading of a disease in a population without immunity, but with disease aware- 
ness. As in the traditional SIS model, infected nodes {I-nodes of fraction [I] of the total node 
number) transmit the disease to an adjacent susceptible node (an S-node of the remaining 
node fraction of size [S]) with rate p, while recovering with rate r. Additionally, S-nodes 
evade infection through selective rewiring (SR) by cutting links to an infected neighbor with 
rate w and rewire it to a randomly selected S-node (double- and self-connections prohibited). 

This basic adaptive network with constant mean degree (k) already displays a rich dy- 
namical behavior in the active phase, featuring small regions of bistable and of stable oscil- 
latory regimes right next to the dominating stable endemic phase lfl9l . All phases, as well as 
the time evolution of both node fractions and the various link densities [AB] (the number of 
AB-lmks, per node connecting nodes of type A and B, A, B £ {S, I}), are approximately 
captured by the pairwise framework proposed in 1191 . While modeling the time evolution of 
the link densities requires tracking triplet densities [ABC] (the number of triplets per node 
with the central node in stage B connected to a A- and C-node (^4, B,C £ {S, I}), the latter 
are approximated by the standard pair moment closure assumption. Then, the evolution of 
the three degrees of freedom of the system in this description is given by 



d[I] 



p[SI] - r[I] 



dt 




(1) 



with (k)/2 = [SS] + [SI] 



+ 111} 



Ifl9l . Since 



(ks) 



2[SS] + [SI] 
[S] 



for the mean degree of the 
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for the I-node subensemble {I- subensemble), Eqs.Q]also offer a coarse-grained view of coe- 
volving network topology. We find that in DE, 



regardless of imposed (k) > 1. Hence in steady state, the S-subensemble will have a larger 
mean degree than the I-nodes if and only if rewiring outperforms infection (see also l43l ). 
This is plausible, since the balance of mean degrees of both subensembles depends on whether 
on average a given Si-link is turned into a SS-link through rewiring or rather into an II-link 
through infection. It follows that (ks) = (kj) in DE if and only if w = p. Remarkably, for 
w = p not only the average degrees of the two subensembles, but also the overall degree dis- 
tributions of the S-nodes and the I-nodes coincide for MC simulations in DE, a feature that 
escapes the simplified description of Eqs. Q] Hence when disease propagation along a SI- 
link and its rewiring are equally probable, a symmetry between the two node subensembles 
ensues. 

The onset of the simple endemic phase, approximated in lfl9l . can be exactly given as 



confirming that rewiring increases the epidemic threshold. Furthermore it can be shown that, 
for any given (k) > 1, w > p + r is a necessary condition for the advent of bistability or 
an oscillatory regime. Consequently a qualitive departure from classic SIS dynamics without 
rewiring is only possible in the pairwise framework if topology change outweighs disease 
dynamics. 

While Eqs. Q] give a coarse-grained description of the adaptive network through various 
averages, the joint steady-state subensemble degree distributions Psj(x, y) of S- and I-nodes 
provide a more detailed probabilistic account and are therefore aimed at by recent models 
0431451 . Since one can extract various network motif densities from the Ps.i(x,y), these 
distributions also encode average lifetimes of any node and link type in a network in DE 
(as do Eqs. [T]by generating triplet densities through the moment closure): In DE, a node 
is on average susceptible for t$ = [S]/ (p[SI]), while the mean lifetime t$i of an STlink 
is determined by considering all of its decay channels, i.e. recovery of the infected node, 
breakup through rewiring and infection of the susceptible partner. Hence 



Similarly, one obtains t$s = [SS]/ (p[SSI]) and tjj ~ 1/ (2r) for the mean lifetimes of 
SS- and II-links, respectively, with all motif densities taken from the network in DE. 

3 Node Cycle 

A new stochastic formalism was proposed in B31 to analytically describe a wide class of 
adaptive networks in DE. Under the assumption that the topology-changing mechanism en- 
sures ergodicity in node-state and node-degree evolution, the node cycle (NC) extracts node- 
ensemble statistics from a single node's long-term behavior in DE, treating its joint degree 
evolution in each state as a random walk. These random walks are coupled by the state- 
change rules and guided by the set /j, of model parameters defining local node dynamics, 
as well as by the set k of correspondence parameters that approximate the influence of 
the network background on local node dynamics in DE. For each /i and k, an analytic ex- 
pression in closed form can be obtained for the transition matrix of the Markov chain in 



(ks) (k z ) = - 



-1, 



w p(k) — r 



w + p + r p ((h) + 1) ' 
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joint degree space associated with a full cycle along the consecutive node states. From its 
Perron-Frobenius eigenvector, whose existence and uniqueness is model-intrinsic, a range of 
probability distributions can be derived describing the stationary state and degree probability 
distributions for arbitrary values of the correspondence parameters. These are for each node 
state A the node's distribution of initial joint degrees $ A , its joint degree distribution P A and 
its stage lifetime distribution T A . Both P* A and T* A are linear tranformations of <P* A in the 
respective node stage A B31 . 

The correspondence parameters must be chosen so as to fulfill any desired global proper- 
ties on the network the NC is unable to account for, such as for instance a fixed global average 
degree. Embedding local node dynamics into a consistent network background imposes addi- 
tional independent constraints on the correspondence parameters. Fulfilling these constraints 
by an appropriate choice of k poses a well-defined, albeit convoluted optimization problem 
(see Sec. [5J. With the optimal choice for k the NC becomes a fully self-contained analytical 
framework for which the long-term dynamics of a single node is made consistent with a given 
network background. Once this has been ensured, the aforementioned distributions represent 
the network's node-ensemble statistics. Important averages can then be obtained from them 
to complement an extensive probabilistic description of the adaptive network in DE: Consid- 
ering lifetime and degree distributions delivers steady-state densities of low-order network 
motifs, as well as mean node-state and link lifetimes H31 . 



4 SIS with Selective Rewiring in the Node Cycle 

For the SR mechanism, the Master equations of the joint degree evolution of a single node 
in a network in DE in stage A (A G {S, I}) can be set up. With [x, y) = Pa(x, y, t\xo, yo) 
being the probability that, having started off with xq susceptible and yo infected neigbors, the 
node under consideration possesses a joint degree of x susceptible and y infected neighbors 
at time t, the node's susceptible S-stage reads as B31 

= { w + r }{( y + l)[x-l,2/ + l] - y[x, y]} 
~py[x,y] +w([x- l,y] - [x,y]) 

+ p s {(x + l)[x + l,y-l}-x[x,y}}, (2) 

and in the infected I-stage as 

^l=r{(y + l)[x-l,y + l}-y[x,y]} 

+ w{(x + 1) [x+ l,y] - x[x,y}} - r[x,y] 

+ p I {(x + l)[x+l,y-l]-x[x,y]}. (3) 

Here the boundary conditions are [—1, y] = [x, — 1] = in both the S- and I-stage, and the 
sets of model and correspondence parameters are /i = {w,p,r} and k = {w,p~s,pi}, respec- 
tively. The correspondence parameters psj approximate the force of infection on susceptible 
neighbors in the two stages, and w the net degree gain through being-rewired-to. The state 
and degree evolution of a single node in a network in DE is then given by a composite ran- 
dom walk described with Eq.|2]in the S- and Eq.[3]in the I-stage. Forpg = w = 0, i.e. with 
frozen network dynamics as described by the NC, the absorbing state of that random walk 
is any [x, 0] (x 6 N) in the S-stage and corresponds to the system's disease-free equilibrium 
also described by Eqs.Q] 

For each [i and k, the steady-state probability distributions Pg j(x,y), <P* S j(x,y) and 
Tg j(t) of a node's degree evolution are calculated from Eqs.|2]and|3] Setting a cutoff for the 
maximum total degree fc max considered in them, this computation is reduced to an eigenvector 
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problem l45l . Distributions Pg j(x, y) yield the mean number (B)a of adjacent B-nodes in 
the A-stage, and similarly one obtains (BC)a, the mean number of (open or closed) triplets 
with the central node in stage A connected to a B- and C-node (A, B,C <G {S, I}). The first 
moments of Pg j(x + y) are consequently given by (k s T ) = (S)gj + (I)s.i an d are the 
mean degrees of the node in the respective stage. The survival functions L* s j(t) of the two 
node subensembles are defined as 

L SJ {t) = l- j T* SJ (t')dt', 
o 

with the mean duration t s 7 of either stage obtained through 

oo 

r*sj = J tT S j(t)dt, 
o 

which for the I-stage yields rf = 1/r as expected B31 . 

Choosing k so that the actual network process in DE is mirrored by Eqs.|2]and|3] distri- 
butions <I>* S j, Pg j, Tg j and L* s 7 describing a node's degree evolution become distributions 
&s,ii Ps,i< Ts,i and Lgj also encapsulating node-ensemble behavior in the respective DE. 
Similarly, NC averages t^j and (kgj) then accurately represent mean subensemble lifetimes 
t g j and mean degrees (kgi)- In DE, the prevalence [I] must be equal to the fraction of the 
time a node typically is infected, and therefore 

ts/t i = (1-[I])/[I] (4) 

must hold, allowing to obtain node-state densities from the NC. 



5 Optimal Correspondence Parameters and Network Stationary 
States 

To tie equilibrium dynamics as described by the NC framework to the actual network process, 
several constraints on the choice of correspondence parameters k need to be imposed. The 
/i and k describing a particular DE can then be used to generate characteristic probability 
distributions within the NC and compare them to the output of individual-based simulations 
(see Sec. 0. Demanding for link conservation (a constant mean degree (k)) yields (k) = 
(Ts(k* s ) + Ti{k*j)) j (t s + ti) and thus 

C „,(l-»±«^) 2 .o. (5) 

For a choice of k that ensures correspondence between NC and network dynamics in DE, 
the three correspondence parameters introduced in Sec. |4]are similar in nature, in that they 
describe the number of infected neighbors of S-nodes, either without additional conditions or 
given that the latter are attached to another particular node type. While then ps.i/p approx- 
imates the mean number of infected neighbors of a susceptible neighbor in the respective 
stage, w/w yields the mean number of infected neighbors of a S-node regardless of the 
neighborhood of the latter. Hence, in the network description, each of these terms can be 
expressed by a mean field composed of low-order motif densities. At exact correspondence, 
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those densities are also given by the NC averages introduced in Sec [4] Consequently, the 
following self-consistency relations should hold regardless of imposed mean degree: 



The NC averages in Eqs. [5] and [6] are convoluted functions of p and k evaluated within the 
NC, and thus for a given adaptive network with fixed p, and (k), each Cj is a cost function in 
K defined by the respective constraint. Because deciding for constant ps.i only approximates 
the interaction of a node's neighborhood with the rest of the network, these four constraints 
may not admit for a solution in k at all, not even for any single constraint. Hence establishing 
correspondence with network dynamics in DE is achieved through identifying optimal k that 
minimize all Ci. 

Due to the different nature of Eq.[5]on the one and Eqs.|6]on the other hand, it is instruc- 
tive to first compute the set of n minimizing Co and of those n minimizing C\ + C2 + C3 
separately. In the absence of degree-status correlations and a characteristic degree imposed 
on the system by minimizing Co, each of the C i = 1,2,3, would be minimal along 
w/w = ps/p = Pi/p — 1- Indeed, MC simulations confirm that the roots of C\ and C2 for 
the motif densities computed from the network in DE verify w/w « ps/p in a wide range of 
p. Thus, for visualization purposes we shall represent the cost functions on the w/w = ps/p 
plane of parameter space (Fig.Q]). 

It can be seen in Fig.QJ) that, the lower w is chosen to weaken degree-status correlations, 
the closer the optimal k for Ci + C2 + C3 is to the diagonal w/w = ps/p = pi/p—1- Since 
candidates for globally optimal correspondence parameters should ensure minimization of all 
cost functions, one has to identify regions of overlap of those k which minimize either only 
Co or Ci to C3. Fixing p and changing the imposed mean degree via Eq.[5]then allows for 
browsing the system's phases as defined by different numbers of overlap regions (Figs.QJ)- 
d)). No region of overlap indicates that the NC cannot be matched to a network in DE, and 
corresponds to parameter values for which the network is in the frozen phase (Fig.QJ)), low 
(k)). One overlap region corresponds to a simple endemic phase (Fig. [TJl), high (k)). Two 
regions of overlap are associated with a bistable phase (Fig.QJ), intermediate (k)): The k 
with larger components represents the stable DE of the stable active branch, whereas the one 
with smaller is associated with the unstable DE of the hysteresis loop leading to bistability 
in Eqs.nHH. 

The approximate values for the k minimizing Co to C3 (i.e the coordinates of overlap 
regions in the w/w = ps/p plane) are used as initial guesses for the minimization of the 
overall cost function Co + C\ + C2 + C3 in full correspondence parameter space. Standard 
minimization routines then yield the optimal correspondence parameters used in the follow- 
ing section to obtain the NC results. 

Since Eqs. [6] also allow for a reverse calculation of optimal k for given set of motif 
densities of a network in DE, they can be used to approximate the expected coordinates of 
optimal k in every overlap region. For stable DEs of a variety of p and imposed (k), these 
predicted coordinates very well match those obtained within the NC. Using the moment 
closure approximation in lfl9l . steady-state densities of triplet motifs can be calculated in 
the pairwise model alongside node-state and link densities. For unstable DEs in the bistable 
phase of Eqs.Q] these densities again give a good estimate for the coordinates of the respective 
overlap region, underlining the validity of the NC approach. 
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a) b) 




w/(w + w) = p s /('p + p s ) w/(w + w) = ps/ ip + ps) 



Figure 1. (Color online) Logarithmic color-coded plots of cost functions C, on w/w = ps/p 
plane with axes rescaled to unit length. Coordinates for low cost function values (< 10~ 5 ) 
are K ensuring good correspondence for the respective constraint, a): C\ + C2 + C3, with 
optimizing k close to, but not on w/w = ps/p- b)-d): Co, with solid lines representing k that 
minimize C\ + C2 + C3 as in a), b): (k) = 3. c): (k) = 5. d): (k) = 7. Model parameters 
w = 0.05, p = 0.008, r = 0.005 and maximum cutoff degree fc max = 50. 



Given the same /i in the NC and MC simulations, the values of imposed mean degree 
(k) that trigger a change in the number of overlap regions coincide with the values of (k) 
observed at corresponding phase transitions in simulations. Furthermore this number is equal 
to the number of stable and unstable DEs predicted by Eqs.Q]for the corresponding phase. 
Hence the NC allows for the localization of all stable and unstable dynamic equilibria for a 
given SIS adaptive network. It is moreover able to structurally explore the dynamics by cor- 
rectly predicting the existence of and the transitions between phases with different numbers 
of DEs. The supercritical Hopf bifurcation in Eqs.[T]that gives rise to a small stable oscillatory 
regime in the bistable phase Ifl9| leaves the number of DEs unchanged. Consequently no in- 
dications for that phase transition were found in the NC, as the framework cannot distinguish 
between the bistable phase's active steady state and its oscillatory regime. Concurrently, the 
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NC is insensitive to the stability of its detected DEs. Yet it strongly hints that unstable DEs 
are physical and not mere model artefacts, and offers, unlike current compartmental models 
and MC simulations, an extensive description of them. 



6 Comparison of Rewiring Mechanisms 

Two modifications of the adaptive SIS model are presented here to showcase the applicabil- 
ity of the NC framework. The first, media-driven rewiring (MR) introduced in P4l . relates 
disease awareness to instantaneous knowledge of the prevalence i through a rewiring rate 
w ■ i (w = constant). Hence in MR, a global time-dependent quantity feeds back to the rate 
of a semi-local rewiring mechanism. The second modification, proposed in [39 ), suggests 
that susceptibles rewire links with a constant rate to a randomly selected node regardless of 
the state of the latter. This blind rewiring (BR) serves as an antipode to the strictly selec- 
tive rewiring put forward in the original model presented in Sec. |2 with parametrizations 
interpolating between those two limiting cases acknowledging partial knowledge of other 
individuals' disease status in a population. 

To model MR in the NC, the additional correspondence parameter i has to be introduced, 
and all rewiring terms in Eqs.[2]and[3]must be rescaled by a factor i. This new correspondence 
parameter is also needed to properly describe an additional mean field necessary for BR. 
Since there a STlink is rewired to an I-node with probability i, both rewiring terms in Eq. [2] 
for the S-stage need to be rescaled by a factor (l — i\ to represent "successful" rewiring to 
another S-node. In Eq.[3]for the I-phase no such rescaling takes place, but its existing rewiring 
term ought to be complemented by adding w (l — i) {[x — 1, y] — [x, y]} to account for the 
"erroneous" linking of S- to I-nodes. With Eqs. [5] and [6] still being a valid set of constraints 
for both MR and BR, i = [I] should hold, resulting in the additional constraint 



through Eq.|U Again the mean lifetimes in Eq.|7]are convoluted functions of model param- 
eters /i = {w,p, r} and correspondence parameters k = {w,p~s,p~i,i} evaluated entirely 
within the NC. Therefore in the two modified rewiring scenarios, cost functions Cq to C4 
ought to be minimized to establish correspondence between the NC and the respective net- 
work model in DE. A simple calculation shows however that in both MR and BR 



holds for every optimal k, so that i is completely determined by the other correspondence 
and model parameters. Therefore the search space for any optimization algorithm aimed at 
fulfilling constraints C to C4 remains three-dimensional. 

Having been provided with the set of model parameters and the mean degree specifying 
an adaptive network, the NC identifies the DEs and the corresponding sets of optimal corre- 
spondence parameters. For all three considered rewiring mechanisms, one can then extract the 
various probability distributions and averages described in Secs.[3]and|4] and compare them 
to the output of MC simulations (Fig. [2j>. Due to recovery being neighborhood-independent 
and happening at a constant rate, Pj (x, y) = &s(x, y) in all three scenarios, i.e the steady- 
state degree distribution of I-nodes is identical to the distribution of their final degrees right 
before recovery B31 . 

For the model parameters and mean degree used in Fig. [2] both SR and BR are in the 
bistable phase fl 191391 . The respective DE is chosen to be in the stable active branch, whereas 
for the same choice of parameters, MR is in its simple endemic phase [44|. There, MR is 
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Figure 2. (Color online) Characteristic distributions of a network in DE for SR (a)-b)), 
MR (c)-d)) and BR (e)-f)). Left column: Degree distributions for S-nodes (circles), I-nodes 
(diamonds), as well as of initial degrees of I-nodes (squares). Right column: Plots of survival 
functions of S-nodes. Solid lines are predictions by the NC. Insets: Comparison of prevalence 
[I]mc taken from MC simulations and [I]nc computed by in the NC (left column), and of 
mean S-lifetimes obtained from MC simulations and the NC (right column), w = 0.05, 
p = 0.008, r = 0.005; (k) = 5, fc max = 80. a)-b): SR (stable active branch of bistable 
phase) with w = 0.095, ps = 0.017, pi = 0.027. c)-d): MR (simple endemic phase) with 
to = 0.12, p s = 0.022, pi = 0.031. e)-f): BR (stable active branch bistable phase) with 
w = 0.17, ps = 0.026, pi = 0.035. MC simulations according to JH with N = 5 • 10 4 
nodes, initial Erdos-Renyi graph and initial prevalence [7]o = 0.6. Statistics were recorded 
at t = 2 ■ 10 4 for 10 4 network realizations. Error bars are smaller than markers. 



equivalent to SR with the rewiring rate rescaled by a factor [I] (the network's steady-state 
prevalence). Consequently degree distributions and survival functions of those two rewiring 
mechanisms resemble each other (Figs.|2^)-d)), with the higher effective rewiring rate in SR 
prolonging the S-stage, lowering overall prevalence and letting both subensembles sample 
large degrees in comparison with MR. The higher the overall steady-state prevalence, the 
more similar are those two rewiring scenarios. 

At very low prevalences however, the SR and BR model essentially describe the same be- 
havior, since they both feature approximately the same rate of successful rewiring to another 
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S-node. If chosen model parameters and mean degree force SR into the bistable phase, there 
will be no low-prevalence active steady state in SR |[T9l . and thus neither in BR. Hence for the 
same set of fj, and (k), the prevalence in a steady-state BR has to be high, and consequently 
successful rewiring happens on a much slower timescale than disease dynamics governed 
by model parameters p and r. Topological separation of node states (through rewiring) is 
consequently outpaced by their randomization (through conventional SIS dynamics), so that 
rewiring effectively randomizes network topology. Much like for very low rewiring rates in 
SR, Poissonian degree distributions around the overall mean degree ensue (Fig. H£)), with 
almost no node-status clustering occuring (as opposed to SR and MR with same \i and (k)). 

The exponential survival function in Fig. |2f) implies exponential lifetime distributions 
Ts(t) of S-nodes for BR, indicating that the force of infection on a S-node, and thus its 
number of infected neighbors, is approximately constant in DE. Of the three processes that 
influence that number, both recovery and successful rewiring have very low rates, the latter 
due to the large prevalence [/]. The rate of the third process - conversion of susceptible 
into infected neighbors - is also small due to the small number of susceptible neighbors 
of a typical node: Neither in its S-stage (due to erroneous rewiring and slow recovery of 
infected neighbors) nor in the preceding I-stage (due to again a modest recovery rate and a 
large number of competing I-nodes at the receiving end of erroneous rewiring) can the node 
agglomorate a large fraction of susceptible neighbors. Consequently the relative change in 
the number of infected neighbors is small, and hence the neighborhood of a S-node exerts an 
almost constant force of infection on it. 

Since the iterative nature of the optimization procedure outlined in Sec.|5]makes the latter 
computionally expensive (see also Sec.[7J>, a relatively small cutoff of total degree £; max in NC 
transition matrices is convenient to identify DEs. However, properly modeling the degree 
evolution of a node with Eqs.|2]and[3]requires setting a sufficiently large fc max . Consequently 
optimal k obtained for low degree cutoffs may provide a slightly blurred correspondence 
between the NC framework and an adaptive network in DE, resulting in small deviations in 
distributions and averages as observed in Fig. [2] 



7 NC Computation and Comparison to other Frameworks 

In the case of large cutoff degrees fc max or a high number n of node states, the procedures 
involved in the NC, as laid out in Sees. [3] to and P31 . can be computationally expensive: 
For each of the n node stages, the Master equation for the respective joint degree evolution 
is set up and solved for arbitrary fi and n. For every k considered, one subsequently com- 
putes n dense stage-transition matrices with ~ fc^" x nontrivial coefficients each. Then the 
unique positive eigenvector of the resulting matrix product yields a probability distribution 
(<P* s (x,y) in the SR scenario) and consequently also steady-state network motif densities 
needed for the cost functions' evaluation at k, allowing for the functions' minimization and 
ultimately identifying optimal correspondence parameters of the DE. A more straightforward 
route to those essential network motif densities is to write down the Master equation for the 
composite random walk through all node stages, i.e. the transition matrix for the full node 
cycle. The resulting equations are equivalent to a stochastic interpretation of the compartmen- 
tal formalism presented in ll43l if time-independent mean fields of link gain in the S-stage 
and of infection acting upon susceptible neighbors are assumed. This transition matrix has 
~ (nfc™ ax ) 2 entries, but is sparse, and its nonzero entries are linear combinations of model 
and correspondence parameters. By arguments similar to those brought forward in B31 . the 
matrix' null space is one-dimensional and yields steady-state degree distributions as well as 
node-state densities (in SR PJ j{x,y) and [I], respectively), directly delivering the crucial 
network densities needed in the quest for optimal k. The computational effort involved in 
evaluating the cost functions in this modified NC thus can be significantly reduced compared 
to implementing the classic NC procedure. The three main modeling frameworks dealing 
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Pairwise Model 


Compartmental Model 


Node Cycle 


transient modeling 


numerical 


numerical 




DE detection 


analytical 


(numerical) 


(analytical) 


DE description 


[A],[AB] 


[A], Pa 


[A], P A , $a, T a 


DE stability 


stable/unstable DEs 


stable DEs detected, 


stable/unstable DEs 




detected, described, 


described 


detected, described 




distinguished 







Table 1. Comparison of frameworks capturing DEs in adaptive networks. Node-state densi- 
ties [A] and subensemble joint degree distributions Pa encode link densities [AB], higher- 
order star motif densities as well as mean node-state and link lifetimes (see Sec.|2]and 1451 ). 
with A and B assuming any value in node-state space. For every node stage A, T| and <P A 
are the stage lifetime distribution and joint distribution of initial degrees, respectively. Iden- 
tifying DEs with the compartmental model is subject to the usual limitations of numerical 
integration, whereas for the NC, cost functions are computed analytically and minimized 
with standard optimization techniques. 



with dynamic equilibria in adaptive networks are contrasted in Tab.Q] To account for correla- 
tions beyond immediate network neighbors, the pairwise model from JT9] can in principle be 
extended to a moment closure at the level of higher-order network motifs, at the expense of 
analytical tractability due to an increasing number of nonlinear ODEs involved. Similarly, the 
accuracy of the compartment model in |43j and the NC can be improved, with computation 
times drastically increasing. 

By construction, the pairwise model is too coarse-grained to distinguish between some 
vital differences in rewiring mechanisms, whereas more local frameworks like the NC or the 
compartmental model can accommodate those changes that do alter the ensuing steady-state 
topology. If for instance a Si-link has been cut in the original SR scenario, adding it between 
two randomly selected susceptibles instead of classically rewiring it to just one implies a 
different random walk in the NC than Eqs.[2]and[3] and hence altered degree distributions 
extracted from them. The equations of the pairwise model however cannot account for this 
change in rewiring mechanism and remain equal to those describing the SR scenario. Like- 
wise, BR reduces like MR to a rescaled SR in Eqs. [TJ whereas in the node cycle and the 
compartmental model, that structural difference between BR and SR is correctly accounted 
for. Moreover the NC in its current form is, unlike the two ODE-based models, not designed 
to capture active phase dynamics other than DEs. But similar to the pairwise model, it can 
indicate a global frozen phase for given model parameters /i: When the NC's different cost 
functions do not display overlapping minima, no self-consistent embedding of single-node 
dynamics in a network in DE exists (see. Fig.QJ))). 

8 Conclusion 

In this paper, we review different frameworks dealing with cyclic node-state dynamics on 
adaptive networks. Focusing on the SIS model with rewiring, we discuss both previous and 
new results that can be derived in the framework of its low dimensional deterministic de- 
scription through the standard pair approximation. By its very nature, this approach cannot 
describe one of the most striking features of the model. This is the fact that for a large param- 
eter regime, simulations on networks show the states of the nodes and the links coevolving 
to produce and maintain a dynamic network topology characterized by well-defined degree 
distributions not only for the global network, but also for the subsets of S- and I-nodes. 

The main focus of the paper is describing these steady-state network statistics, for they 
are essential to understanding the relation between dynamics and network structure and to 
apply these ideas to real world examples. We briefly present the ideas of a method that has 
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been proposed recently to obtain an analytic description of the dynamic equilibria of adaptive 
networks. We then apply this method, called the node cycle, to the SIS model with three 
different types of rewiring. These three different rewiring schemes cover a whole range of 
assumptions about the less well-established "social" dynamics of the model. 

A key point of the node cycle method is the determination of the "correspondence pa- 
rameters", whose role is to embed the general node-cycle description into the dynamics of a 
particular network through the fulfilment of self-consistency conditions that express the net- 
work's topology in dynamic equilibrium. We discuss in detail how to determine the optimal 
values for these correspondence parameters within the node cycle framework. For a given set 
of model parameters, this self-consistency approach detects and describes all dynamic equi- 
libria predicted by the pairwise model, notably including a second equilibrium in the bistable 
phase that can be identified with an unstable active steady state of the system. 

The node cycle method allows the determination not only of the degree distributions for 
all node stages, but also of their distributions in initial degrees, their lifetime profiles and of 
the densities of low-order network motifs in the steady state. We contrast its scope of de- 
scription of dynamic equilibria with that of both the compartmental and the pairwise model 
and identify areas of competence for each of the descriptions. The computational effort in- 
volved in the node nycle model is laid out, and an alternative, faster route to aforementioned 
steady-state characteristics is given. For the three rewiring schemes that we considered, the 
results for the steady-state degree distributions and lifetime profiles derived from the node- 
cycle method are in very good agreement with the results of simulations on networks. As 
to the comparison of the different rewiring schemes, we find that modifications of the stan- 
dard rewiring rules with very similar functional forms can have either negligible or dramatic 
effects on the resulting network topology. 

In conclusion, we review several frameworks dealing with dynamic equilibria in adap- 
tive networks, elaborate on the pairwise formalism of Gross et al. that models SIS dynamics 
on adaptive networks, apply the node-cycle method and show that it can be used to iden- 
tify steady states. We subsequently describe the latter through degree distributions and other 
steady-state statistics. In future work we intend to explore the method further by applying it 
to other processes on adaptive networks beyond a binary state space. 
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